Aharonov-Bohm oscillations in the local density of states 
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The scattering of electrons with inhomogeneities produces modulations in the local density of 
states of a metal. We show that electron interference contributions to these modulations are affected 
by the magnetic field via the Aharonov-Bohm effect. This can be exploited in a simple STM setup 
that serves as an Aharonov-Bohm interferometer at the nanometer scale. 



Scanning tunneling microscopy (STM) provides a pow- 
erful tool to measure the density of states py. The 
scattering of electrons with inhomogeneities is known to 
modify the corresponding density of states producing lo- 
cal modulations. These modulations were first probed 
by STM with atomic resolution in [2 , where standing- 
wave patterns in the local density of states (LDOS) aris- 
ing from scattering of electrons with impurities and step 
edges were observed on the surface of Cu(lll). These 
patterns encode information about the electronic prop- 
erties of the corresponding system. For example, in the 
case of a two-dimensional electronic system, a point de- 
fect produces oscillations in the LDOS with wavevector 
2kF for very low bias volatge [2 . The period of these os- 
cillations changes as a function of the bias voltage, from 
which one can infer the spectrum of the electrons. In 
the case of a superconductor, there can be additional 
features reflecting the symmetry of the superconducting 
gap [3 . Moreover, atomic manipulation permits the en- 
gineering of surfaces, making it possible to confine elec- 
trons into the so-called quantum corrals [4 . This has 
attracted great attention since it permits the study of 
lifetime effects, Kondo physics, single-atom gating, etc. 
[5 . Very recently, open nanostructures have been shown 
to be amenable for quantum holographic encoding [6] . In 
the following we show that under suitable conditions, the 
LDOS exhibits oscillations due to the magnetic field that 
can be interpreted as due to the Aharonov-Bohm effect. 
Therefore the STM setup can be designed to serve as an 
Aharonov-Bohm interferometer at the nanometer scale. 

For the sake of concreteness we consider the close- 
packed surface of a noble metal where the so-called 
Shockley surface states form a two-dimensional nearly 
free electron gas. Two atoms deposited on top of this sur- 
face can be modeled as two point scattering potentials for 
the surface electrons. This forms the simplest Aharonov- 
Bohm interferometer, where the role of the two different 
paths in conventional Aharonov-Bohm setups is played 
by the two scattering paths shown in Fig. [l] In the pres- 
ence of a magnetic field, electrons scattering along these 
loops pick up different phases depending on whether the 
scattering is clockwise or anticlockwise (the two paths be- 
ing connected by time-reversal symmetry at zero field). 




FIG. 1: STM interferometer, r represents the position of 
the STM tip on the surface and ri and r2 two impurities. 
The LDOS measured by the STM tip contains interference 
contributions due to electrons traveling along the two paths 
shown in the figure. The magnetic field affects this interfer- 
ence via the Aharonov-Bohm effect, producing oscillations in 
the LDOS that is measured by the tip. 

This affects the interference contribution to the LDOS 
measured by the STM tip, which eventually exhibits os- 
cillations as a function of the magnetic flux that passes 
through area enclosed by the above paths. The physics 
of these LDOS oscillations is similar to the effect of the 
magnetic field on weak localization in disordered two di- 
mensional metals [7 . 

The (i/(r, V)/dV maps obtained experimentally are de- 
termined by the LDOS of the sample A/'(r,a; = eV)^ the 
tip density of states, and the tunneling matrix elements 
^ ^ . In the Tersoff-Hamann approximation with a con- 
stant tip density of states dl/dV is proportional to N 
[S'.'S^. However, the oscillations of the LDOS that we dis- 
cuss are picked up by the experimental dl/dV even if the 
above proportionality is lost due to voltage dependence 
of the tunneling matrix elements or due to variation of 
the tip density of states with energy. The LDOS of the 
sample can be obtained from the corresponding retarded 
Green's function as 

N{t,u;) = --lmG^{T,T;u) (1) 
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(the factor 2 is due to spin degeneracy). For a two- 
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dimensional free electron gas the Green's function is 



The presence of impurities can be modeled by a term 



(2) 



Here N = m/{27TJi?) is the density of states of the elec- 
tron gas per spin, Hq'^ is the zeroth-order Hankel func- 
tion of the first kind and k{Lj) is given by the dispersion 
relation k{u) = kpil + Cc;//i)-^/^, where /i = h^k\/ {2m). 
For large distances [r ^ k~^{(jo)] we have 
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in the Hamiltonian of the system, where U is the scatter- 
ing potential associated with the impurities and p is the 
electronic density. Following a perturbative approach [5j 
the Green's function can be expressed as G = Go + (^G, 
where Go is the Green's function in the absence of impu- 
rities and 



5G{v, r') = j dv"Go{v - v")U{v")Go{v" - r') 

dY"dv"'Go{v - y")U{y")Go{v" - r"')C/(r"')Go(r"' - r') + . . • 
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The dependence on the frequency is dropped since the 
scattering is assumed to be elastic. The latter quantity in 
([5| contains the interference contributions to the LDOS 
we are interested in (i.e. from terms second order and 
higher in the scattering potential). 

Two identical point impurities are described by the 
scattering potential 



C/(r) = C/o[(5(r-ri) + 5(r-r2)]. 



(6) 



In this case the quantity (^G(r, r) naturally contains two 
types of terms. On one hand, there are (additive) terms 
in which the scattering with the impurities is produced 
separately. The Fourier transform of this contribution 
has quasiparticle peaks which contains information about 
the dispersion of the electrons (see e.g. [3]). This con- 
tribution, however, plays no role in the physics that we 
intend to study. On the other hand, there are terms 
involving scattering with both the impurities. Among 
the latter terms, there are processes in which the semi- 
classical scattering paths enclose a finite area as in Fig.[T] 
It is important to note that such loops occur in pairs 
connected by time reversal symmetry (i.e., clockwise and 
anti-clockwise). For example, to the lowest order in the 
impurity potential, the contribution to the Green's func- 
tion due to closed loops reads 
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In the absence of a magnetic field the two terms in this ex- 
pression are equivalent (due to time reversal symmetry), 
and therefore give the same contribution to the LDOS 
(it can be said that the interference is constructive). In 



the presence of a magnetic field, however, time reversal 
symmetry is broken. Then electrons travelling clockwise 
and anti-clockwise along the above loop acquire different 
phases, so the subsequent interference is affected. This 
results in Aharonov-Bohm oscillations of the LDOS as 
we shall see explicitly below. 

Next we demonstrate that the simple geometrical pic- 
ture above remains unchanged when higher order terms 
in the impurity potential are taken into account. At 
higher order we have to deal with the following additional 
ingredients: (i) multiple scattering at the impurities and 
(ii) multiple scattering where semi-classically the parti- 
cle goes back and forth between the two impurities. The 
first type of multiple scattering can be easily taken into 
account by replacing each of the scattering potentials by 
their respective T-matrices. That is [9^, 
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As regards the second point, we note that 5Gioop(r, r) 
is entirely due to processes where the path between the 
impurities is traversed an odd number of times (otherwise 
the scattering path does not enclose a finite area). In 
the second order contribution ([7|), for example, the path 
between the impurities is traversed once. At 0([/o ) there 
are contributions where the path is traversed three times, 
and so on. Taking all this into account, we obtain 



.5Gioop(r,r) = W^Go{v - ri)Go(ri - r2)Go(r2 - r) 

+ (1^2), (9) 
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It is worth noticing that in ah the processes that finahy 
give rise to ([9| we are actually dealing with the same 
area, since to go back and forth along the same line, for 
example, does not change the area of the resulting loop. 
This is the simple reason why the non-trivial phase rela- 
tion between the clockwise and the anticlockwise paths 
in the presence of a magnetic field is not washed out by 
multiple impurity scattering. 

Let us now consider explicitly the influence of the mag- 
netic field. In the low-field regime (see below) we can use 
the semi-classical approximation for the electron Green's 
function 

Go(r - r') = exp A(l) • Goo(r - r'). (11) 

Here Goo(r — represents the Green's function in the 
absence of magnetic field, A is the vector potential 
(B = V X A), <l>o = h/(2e) is the fiux quantum, and 
the integral is along the straight line connecting r and 
r^ The magnetic field then enters the interference con- 
tributions to the LDOS via complex factors e^*^*^/*^°, 
where ^ is the magnetic fiux through the area enclosed 
by the corresponding scattering path (the different signs 
of the phase corresponds to anti-clockwise and clockwise 
line integrals respectively). As a result, the LDOS can 
be written as 

N{v,uo) = 7VB=o(r,cc;) +7Vioop(r,cc;)[cos(7r<I>/<I>o) - 1], 

(12) 

where Nb=o is the zero-field total LDOS and 



A^ioop(r,cj) = --ImJGioop(r,r;cj) (13) 
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represents the (constructive) interference contribution 
due to all closed paths in the absence of a magnetic field 
[computed from (|9| with Goo]- This interference pro- 
cess picks up a non-trivial phase in the presence of a 
magnetic field, giving rise to oscillations in the LDOS as 



in (12), which is a manifestation of the Aharonov-Bohm 
phenomenon [10 . This can be revealed by varying ei- 
ther the magnetic field or the relative position between 
the STM tip and the impurities since ^ changes in both 
cases. We note that, irrespective of the position of the 
STM tip, the magnitude of the correction to the LDOS 
reduces with B for low magnetic fields. As in the case 
of negative magnetoresistance in weak localization, this 
is due to the fact that magnetic field induces destructive 
interference between the contribution of the two semi- 
classical paths. 



Next we discuss the limitations of our calculations 
and the feasibility of our proposal to use the STM as 
an Aharonov-Bohm interferometer at a nanometer scale. 



The expression (12) for the LDOS has been derived 



within a semi-classical approach, so it holds as long as 



Eq. (11) can be used for describing the infiuence of the 



magnetic field on the electron system. This is possible if 
the magnetic field is such that the Fermi wavelength is 
much smaller than the Landau orbits: 
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The magnetic field needed to observe a complete 
Aharonov-Bohm oscillation can be estimated as 
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where d is the characteristic distance in the setup (i.e. 
the distance between the impurities and/or between the 
impurities and the position of the tip. For d ~ 40 — 20 nm 
this field is B27r ~ 5 — 20T, which corresponds to Lan- 
dau orbits ~ 11 — 6nm. For Cu(lll) and Ag(lll) 
Xf = 2.95 nm and 7.6 nm respectively, and therefore the 



condition in Eq. ( 14) for the semi-classical approximation 



is valid. Furthermore, inspite of the fact that the inter- 
ference signal is long-range in the sense that TVioop decays 
as l/d in a power-law fashion [see (|3])], in reality there are 
dephasing processes that introduce an extra attenuation 
(thermal dephasing for example), and which, for simplic- 
ity, have not been taken into account in the current com- 
putation. Experimentally, in fact, the impurity-induced 
variations in the LDOS are typically observed up to dis- 
tances of the order of a few times the Fermi wavelength, 
say IOXf [II]. Therefore the characteristic distance in 
our setup must be d < IOXf- For d ^ 40 — 20nm we then 
need A^? > 4 — 2nm, which still can be smaller than at 
5 — 20 T. Thus we find that the semiclassical interpreta- 
tion is good to describe the first few periods of the LDOS 
oscillations. The first corrections to our results will be 
due to the curvature of the classical trajectories, which 
can still be described within the semiclassical picture [12] . 
To describe sufficiently high periods, however, one has to 



go beyond ( 11 ) and consider the inffuence of the magnetic 



field within a Landau-level approach. These periods im- 
ply magnetic fields considerably high (> 20 T), so we do 
not develop this latter approach here. 

The Aharonov-Bohm physics reveals also in the spatial 
variations of the LDOS for a fixed magnetic field. In 
our setup ^ varies only in the direction perpendicular to 
the line connecting the two impurities. In consequence, 
the STM scans along this direction will show a periodic 



envelope due to the cosine factor in Eq. (12). With two 
impurities separated 20 nm, fields of 5 — 20 T give rise to 
periods of 80 and 20 nm respectively for such an envelope 
(see Fig. [2|. 

It is also worth mentioning that the Zeeman splitting, 
neglected so far, has a trivial infiuence on the Aharonov- 
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FIG. 2: Expected STM patterns for two impurities 20 nm apart on the Ag(lll) surface after subtraction of the B = signal. 
The horizontal stripes are produced by the Aharonov-Bohm effect in the LDOS, whereas the remaining elliptic features are due 
to the interference contributions described by Moop [with Af = 7.6 nm and real VK-matrices as defined in Eq. (10)]. 



Bohm oscillations if the spin is conserved in the scat- 
tering process. Moop in Eq. (12) actually results from 



contributions associated with the two spin polarizations. 

It can be written as TVioop = U^lop + ^/oop) 

spin is conserved. The eventual difference between N^^^^ 

and A^ioop due to the Zeeman splitting can be probed 
by means of spin polarized STM. However this difference 
does not alter the oscillatory behavior of the LDOS de- 
scribed above. The situation is more subtle if the spin 
can flip during the scattering process as a result of e.g. 
spin-orbit coupling. This can affect multiple scattering 
in a nontrivial way [13 , and can give rise to the analog 
of the anti-localization phenomenon if spin components 
are measured separately. 

In summary, we have shown that magnetic field affects 
electron interference contributions to the local density of 
states via the Aharonov-Bohm effect. This can be ex- 
ploited in building STM devices that serve as Aharonov- 
Bohm interferometers at the nanometer scale. We have 
illustrated this possibility for the close-packed surface of a 
noble metal with two atoms adsorbed on top. The role of 
these atoms consists in creating strong enough scattering 
potentials, which can also be produced, for example, us- 
ing additional STM tips. The implementation of this new 
functionality into the STM technique might broaden its 
applications notably, offering new perspectives for STM 
studies of the fundamental properties of surfaces and un- 
derlying systems. 
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